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We present an improved calculation on the pionic twist-3 distribution amplitudes (j>p and < 

By adding all the uncertainties in quadrature 
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which are studied within the QCD sum rules. 

is found that (g) = 0.248™, (t P ) = 0.262™, (&) = 0-102™ and (£) 



n 094+ 028 
u.uyi_ 020 . 



Furthermore, with the help of these moments, we construct a model for the twist-3 wave functions 
■tpp i<T (x,li±), which have better end-point behavior and are helpful for perturbative QCD approach. 
The obtained twist-3 distribution amplitudes are adopted to calculate the B — >• ir transition form 
factor f£ within the QCD light-cone sum rules up to next-to-leading order. By suitable choice of 
the parameters, we obtain a consistent fi with those obtained in the literature. 



o 

(N 
C 

>—> 

On 



43 ■ 

Or 
i ' 

<d : 
43 



> 

OS 

in 

m 



x 



PACS numbers: 11.55.Hx, 14.40.Aq, 12.38.Aw 



I. INTRODUCTION 

Distribution amplitude (DA) shows the momentum 
fraction distributions of partons in hadron in a particu- 
lar Fock state, which posses an important component for 
QCD factorization theory [lj. It is convenient to arrange 
DA by its different twist structures such that when the 
energy scale of the process is large enough, one can safely 
neglect those higher power suppressed contributions from 
higher twists. Due to its simpler structure, the leading- 
twist DA has attracted much attention in the literature 
1, 2]. Very recently, the new BABAR data on the pion- 
photon transition form factor [3j arouses people's new 
interests on the pion twist- 2 DA [J]. 

However the higher-twist DAs are much more involved 
0-0 > which describes either contributions of the trans- 
verse motion of quarks (antiquarks) in the leading-twist 
components or contributions of higher Fock states with 
additional gluons and/or quark-antiquark pairs or etc.. 
In certain cases, higher twist structures especially the 
twist-3 DAs may provide sizable contributions and there 
are usually the main uncertainties for the theoretical es- 
timation. For example, as for the case of pion electro- 
magnetic form factor and B — > n transition form factor, 
most calculations give large twist-3 contributions that are 
even dominant over that of the leading twist in the inter- 
mediate or even large Q 2 region. It should be pointed out 
that those large contributions of twist-3 DAs are usually 
based on the asymptotic behavior of (fE, i.e. (j)p S = 1. 
And with such a naive asymptotic behavior of twist-3 
DA, the end-point singularity can not be effectively sup- 
pressed.which inversely leads to a large twist-3 contri- 
bution |a-(l2l|. While, by taking the kx factorization ap- 
proach [131 ] and by constructing a twist-3 wavefunction 
model based on the DA moments derived by Refs. [l4l . ll5j . 
Refs.flrl [l7j show that even though in the intermediate 
energy region, both leading twist and twist-3 DAs can 



provide sizable contributions, however the contributions 
from twist-3 ones are really power suppressed for large 
Q 2 region, and then the normal power counting rule is 
retrieved. So a twist-3 DA with a better end-point behav- 
ior other than the asymptotic one shall lead to a better 
understanding of these form factors. 

At the present, the pionic twist-3 structures are far 
from affirmation 0, El El EE El , and it would be 
interesting to do further studies on the twsit-3 DAs / 
wavefunctions. In the present paper, we shall first make a 
study on the pionic twist-3 DA moments within the QCD 
sum rules, and then make a discussion on its physical 
effects by constructing a reasonable wavefunction model 
and by taking B — > 7r transition form factor as an explicit 
example. 

The remaining parts of the paper is organized as fol- 
lows. In Sec. II, we present the calculation technology 
for the twist-3 DAs tfSE and within the framework of 
QCD sum rules. In Sec. Ill, we present our numerical 
results for (ft™ and </>J, where the DA moments together 
with their uncertainties are discussed in detail. Based on 
the derived DA moments, we construct the wavefunction 
models for ip* and tp™ in Sec. IV, and then make a discus- 
sion on their rationality by further applying them to deal 
with the B — > tt transition form factor within the QCD 
light-cone sum rules up to next-to-leading order (NLO). 
The final section is reserved for a summary. 



II. CALCULATION TECHNOLOGY FOR THE 
TWIST-3 DAS 4>l AND 4>l 

The pionic twist-3 DAs and </>J are defined as El 
(0 \u{x)i l5 d{-x)\ 7T-) - I due^ x Ul(u) (1) 



and 
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due i( ^ x \q^x v - q v x^)<t>l{u), 



(2) 



where £ = (2u — 1), and due to the equations of mo- 
tion of the quarks inside pion [l9| , = — m + md i whose 
value is 1.6 ± 0.2GeV at the scale 1 Gev'Q. Normally, 
these two DAs can be expanded by the conventional 
Gegenbauler expansions, = J2n (£p n ) ^2niCl an<1 

m) = 3(1 - 0/4E„ <e) C^ 3 (0, where C^ a ' 3/a (0 



are Gegenbauler polynomials and (C 2 ,^) are the so-called 
Gegenbauler moments. As for the pionic case, because 
of the chiral symmetry, only even moments' terms are 
non-zero. 

The Gegenbauler moments can be calculated under the 
QCD sum rules. For such purpose, one can define the two 
correlation functions as 



(z ■ q) 2n I^°\q 2 ) = ~i d 4 xe lq - x (O T{d(x) lb (iz ■ D) 2n u(x), u(0) 75 d(0)} 



(3) 



and 



- i(q„z v - q v z,){z ■ <7) 2n 4 2n '°V) = -i fd 4 xe l " x (o \T{d(x)a„ v75 (iz ■ D) 2n+1 u(x), u(0) l5 d(0)}\ o) . (4) 



Following the same calculation technology as described in 
detail in Refs.[l4l. [l5j. i.e. the QCD sum rules under the 



background field approach [2CH23l |. we obtain two sum 
rules for the moments of and with the condensates 
up to dimension six, i.e. 
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M 2n + 1 
2n+l(^G 2 ) 



1 



in 2 2n + 3 



M 2 



24 



M 4 



167T .„ 9 „. ( 

(8?i 2 - 2p 

81 1 M 6 



- V 

qq) 



2n + 1 (m u + m d ) {qq) 
2 M 4 



(6) 



where M is the Borel parameter, and are contin- 
uum threshold, which are usually taken to be around the 
mass square of the first exciting state 7r'(1300) of pion. 
For clarity, we will take s^ a = 1.59, 1.69 and 1.79 GeV 2 
to do our discussion. For the non-perturbative vacuum 
condensate, we take [H|: (uu) = (dd) ~ -(240 MeV) 3 

and ^— G 2 ^ = 0.012 GeV 4 , which are at the renormal- 

ization scale 1 GeV and shall be evaluated up to the 
concerned scale by using the leading order anomalous 
dimensions of the condensates. The pion decay con- 
stant, the pion mass and the a s are taken as the cen- 
ter values of Refs.[H H|, i.e., U = 0.1304 GeV and 



m n = 0.1396 GeV and a s (M z ) = 0.1184. The renormal- 
ization scale for the present case is taken as fi = M . 

As a cross check, we make a comparison with the sum 
rules derived in Refs.[l4l.ll5j. Firstly, our sum rules agree 
with the corresponding sum rules derived by Ref. (l5| , ex- 
cept for the two coefficients of the dimension- four (yG 2 ) 

and dimension-six matrix element ( y/o^qq) 2 . It is found 
that such a difference is caused by the two Feynman 
diagrams shown by FIG.([T|), which are missed by the 
authors of Re f. 11511 . According to the background field 
approach [20-23], the Feynman rules are derived un- 
der the Schwinger Gauge or 'the fixed-point gauge' [25| . 
And because of the gluon's or quark's equation of mo- 
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FIG. 1. Feynman diagrams that are missed in Ref.ll5j|. where 
the gluon and the quark condensates are depicted as crosses, 



and the big dot stands for the currents 
and d(x)a lll yj 5 (iz ■ D) 2n+1 u(x) for q 



d(x)j s (iz ■ D) 2n u(x) 
] and 4>& respectively. 



TABLE I. Uncertainties of the second moment , where 

/^(lGeV) = 1.4, 1.6, 1.8 GeV respectively. 



<(GeV 2 ) 


1.59 


1.69 


1.79 


M 2 (GeV 2 ) 

(0(1.4) 
(0(1.6) 
(O(i.s) 


[0.696, 0.874] 
0.315 ±0.002 
0.242 ± 0.002 
0.191 ±0.001 


[0.696, 0.905] 
0.324 ± 0.005 

0.248 ± 0.004 

0.196 ±0.003 


[0.696, 0.937] 
0.334 ± 0.009 

0.255 ± 0.007 

0.202 ± 0.005 



tion, only those interactions that have at most one back- 
ground field coupling can have contribution to the inter- 
action vertex. As for the two Feynman diagrams shown 
by FIG.([T]), the two background fields are attached to 
a big dot, which stands for the two complex currents 



d(x)j 5 (iz- D) 2n u{x) 
for d>Z and 



and d{x)<7^ v "fs(iz ■ D) 2n+1 u(x) 

respectively. Since the two background 
fields can attach to different places of the two complex 
currents, these two Feynman should be taken into con- 
sideration. Secondly, Ref. [l4[ only presents the sum rules 
for (pp. And, one may find a typo error in Ref. [l4j . i.e. 
the term before the dimension-three quark condensate 
< qq > and dimension-six quark condensate < qq > 2 
should be " ± " other than " — " . 



III. NUMERICAL RESULTS FOR THE 
MOMENTS OF <j>* AND 

Basing on the QCD sum rules (JSJ and dU), we discuss 
the Gegenbauler moments of tp™ and </>^ . For the purpose, 
we adopt the usual criteria to fix the Borel window for 
the first two moments (Cp^)j i- e - the continuum contri- 
bution is less than 30% of the total dispersion integration 
and the dimension-six condensate contribution is no more 
than 30%. It is found that the continue contribution in- 
creases with the increment of M 2 and the contribution 
from the dimension-six condensate term increases with 
the decrement of M 2 , then possible Borel windows can 
be obtained as required. 

All the obtained Borel windows and the corresponding 
second and fourth moments of and (/>J are collected 
in TABs.(lIininilIV|) respectively. TABs.(|]jnjnijIV|) in- 



TABLE II. Uncertainties of the fourth moment (£p} , where 
MlGeV) = 1.4, 1.6, 1.8 GeV respectively. 



<(GeV 2 ) 


1.59 


1.69 


1.79 


A/ 2 (GeV 2 ) 

(0(1.4) 
(0(1.6) 
(0(1.8) 


[0.958, 1.052] 
0.335 ± 0.004 
0.257 ± 0.003 
0.203 ± 0.002 


[0.958, 1.080] 
0.342 ± 0.004 
0.262 ± 0.003 
0.207 ± 0.002 


[0.958, 1.108] 
0.349 ± 0.004 
0.267 ± 0.003 
0.211 ±0.002 



TABLE III. Uncertainties of the second moment (Cct) 
where /^(lGeV) = 1.4, 1.6, 1.8 GeV respectively. 



<(GeV 2 ) 


1.59 


1.69 


1.79 


M 2 (GeV 2 ) 

(W(1.4) 
(W(1.6) 
(W(1.8) 


[0.439, 0.798] 
0.129 ±0.013 

0.099 ±0.010 

0.078 ± 0.008 


[0.439, 0.832] 
0.134 ±0.017 

0.102 ±0.013 

0.081 ±0.010 


[0.439, 0.867] 
0.139 ±0.021 

0.106 ±0.016 

0.084 ±0.013 



dicate that 

• The second and fourth moments of <j)p are slightly 
affected by the Borel parameter M 2 . This shows 
that the present sum rules ([5]) is reasonable. More 
explicitly, it can be found that the second and forth 
moments of <^p, the uncertainty caused by M 2 is 
less than 2%. For sum rules ©, it is noted that 
for the second moment of 0£ , the uncertain caused 
by M 2 reaches up to a somewhat larger value ~ 
13%. While for higher moments, it is found that the 
uncertainty caused by M 2 is back to be less than 
2%. If we require the flatness of the moments versus 
M 2 as an extra condition to constrain the Borel 
window then the large uncertainty for the second 
moments can be reduced. Here to make a consistent 
analysis over all the moments of (f>p and <f>* , we do 
not imply any further constraint for To show 
these points more clearly, we draw the curves of 
the second and forth moments of a versus M 2 in 
FIGs. (I2I3I) . which shows that the moments (except 
for the second moments of <jj^) are almost fiat in 
the allowable Borel window. 



TABLE IV. Uncertainties of the fourth moment 
where /x^(lGeV) = 1.4, 1.6, 1.8 GeV respectively. 



s£(GeV 2 ) 


1.59 


1.69 


1.79 


M 2 (GeV 2 ) 

(^ CT )(1.4) 
(^ CT )(1.6) 
(^ CT )(1.8) 


[0.807, 0.988] 
0.120 ±0.001 

0.092 ±0.001 

0.072 ±0.001 


[0.807, 1.017] 
0.122 ±0.001 
0.094 ± 0.001 
0.074 ± 0.001 


[0.807, 1.047] 
0.125 ±0.001 

0.096 ± 0.001 

0.076 ± 0.001 
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FIG. 2. The second and forth moments of twist-3 DA cj>p ver- 
sus the Borel parameter M 2 corresponding to fi^ = 1.6GeV. 
The solid line, dashed line and dotted line correspond = 
1.59, 1.69, 1.79GeV 2 , respectively. 



FIG. 3. The second and forth moments of twist-3 DAs cj)Z ver- 
sus the Borel parameter M 2 corresponding to fj,^ = 1.6GeV. 
The solid line, dashed line and dotted line correspond = 
1.59, 1.69, 1.79GeV 2 , respectively. 



• The Borel windows will be broadened with the in- 
crement of the continue threshold s^' a . And ac- 
cordingly, the second and the fourth moments of 
0p CT shall increase with increment of s^ a . When 
increasing s p :' 7 by 0.10 GeV, the second and fourth 
moments of <$L shall be increased by 3% — 5%. 



where the center value are obtained by taking all 
the input parameters to be their center values, i.e. 
S P' a = 1.69 GeV 2 and ^{lGeV) = 1.6 GeV. And 
the uncertain of a particular parameter is obtained 
by fixing other parameters to be their center ac- 
cordingly. 



The twist-3 DA moments are greatly affected by the 
value of /in, which provides the dominant uncer- 
tainty for the present sum rules. So a better deter- 
mination of shall greatly improve our knowledge 
on twist-3 DA. All the twist-3 DA moments are de- 
creased with the increment of /v- When increasing 
/Itt by 0.2 GeV, the second and fourth moments of 
<t>* shall be decreased by 20% - 25%. 



As a summary, by adding the uncertainties in 
quadrature, we obtain = 0.248±°;°^, = 



2 62+ - 080 
u.zoz_ 055 , 

n r)94+ - 028 



(g) = 0.102 



+0.035 
-0.025 



-0.052' 

and (£) 



IV. MODELS FOR Vp AND AND THEIR 
APPLICATION FOR B -> tt FORM FACTOR 



• By taking all uncertainty sources into considera- 
tion, we obtain 



/f 2 \ = 
'Pi 



and 



0.248 ± 0.004| M 2 
0.262 ±0.003 



+0.007 1 
-0.006 Is* 

+0.005 I 
M 2 -0.005 I si 



+0.076 I 
-0.052 IMtt ' 

+0.080I 
-0.055 IMttJ 



(£) =0.102 ±0.013| M2 ±° ™i\ sz t°;°oliV (7) 



0.094 ±0.001 



M 2 



f0.002i +0.028 1 
-0.002lsJ -0.020lMn 



(8) 



A. Models for and i/C 

The wavefunction is the key component for the kx fac- 
torization approach [13j , which is helpful to suppress the 
end-point singularity in certain processes. And it would 
be useful to know the properties of the twist-3 wavefunc- 
tions ipp and 

The wavefunction is purely nonperturbative, so it can 
hardly be determined from the first principle of QCD. 
Generally, the pion wavefunction and its correspond- 
ing DA is related through the relation, <j)l a {x, ///) = 



5 



J\k x \<nf T^^p^^ x ' k± -^ where fj, f is the factorization 
scale that is around O(lGeV). So, on the other hand, 
if we have known the DA moments well, then we can 
inversely obtain some valuable properties of the wave- 
function. 

Our model for the pionic twsit-3 wavefunctions is based 
on the assumptions that I) its longitudinal distributions 
is determined by its DA moments as derived by the sum 
rules ([5]) and ©. And its longitudinal behavior is dom- 
inated by the first two Gegenbauler moments; II) its 
transverse momentum dependence is determined by the 
so-called BHL-prescription [2^, [3(| , which is obtained by 
using the connection between the equal-time wavefunc- 



tion V'c.m.(q_L) in the rest frame and the light-cone wave- 
function ipLc{ x i k_i_) i n the infinite momentum frame, i.e. 

^c.m.(qx) <->• fee ( 4x(i-x) ~ ™ 2 )' wnere m stands for 
the light constitute quark mass. It is found that the 
the transverse momentum dependence is just in the ex- 
ponential form of the off-shell energy of the constitute 
quarks, which agree with the Brodsky and Teramond's 
holographic model that is obtained by using the anti-de 
Sitter / conformal field theory correspondence [3l|, |32| ■ 
And, the pionic twsit-3 wavefunctions ipp(x,\i±) and 
-0J(x,k_L) take the following form 



-0p(x,kj 



1 + B v Cl ,2 {2x - 1) + C p Cl /2 (2x - 1) 



and 



W(x, k x ) = [l + B a Cl /2 {2x - 1) + C a Cl' 2 {2x - 1) 



x{\ 



-exp 



ki 



x{l 



-exp 



Sf3 2 p x{l - a:) 



/ m 2 + kj 



(9) 



(10) 



The parameters A Pt<T , B p<Tl C PiCr and /3 PjCr can be 
determined by the average value of the transverse 
momentum (kj_) , the wavefunction normalization 

Jo dx i|k x |< w Tsfr^pA^ k -L) = !> and the nrst tw0 DA 
moments determined by the last section. The average 



value of the transverse momentum square is defined as 



/ dxd 2 'k± 


k 2 


b^.aO^k 


Ol 2 


J dxd 2 k 


± \% 





and (kl)^ 2 = 0.350GeV ||. Using the relation be- 
tween the DA and the wavefunction, we obtain 



and 



B a Cl /2 {2x 



1) 



' C^cJ (2x 



l) + C a Cl /2 {2x-l) 



^ 6XP 



8^x{l-x) 



2ir 2 



-exp 



8/3 2 x(l - x) 



1 — exp 



1 — exp 



SP 2 x(l - x) 



8/3 2 x{l - x) 



(12) 



(13) 



The Gegenbauer moments (£, p n a ) for the DAs defined 
by Eqs. (|12ll3l) at the scale fif can be defined as 

/t2n\ I ! 1 Q dxr p {x^ f )C 1 2 j 2 (2x-l) 



and 



tidx[Cl' 2 (2x-\W 



U2n\ | Jo dx(j>Z(x,fi f )C^ 2 (2x- 1) 

tidx6x(l-x)iCtL 2 (2x-lW [ ] 

With the help of the above formulae, we can determine 



the wavefunction parameters for varying DA moments, 
which are presented in TAB |V] and TAB ED 

We present 0£ and <f a in FIGs. (jiI5|) . Our DAs, i.e. 
Eqs. (|12ll3p . are drawn by solid lines and are derived by 
setting the second and fourth moments to be their cen- 
ter values shown in TABs. (|IIIIIIIIIIV| ). As a comparison, 
we also present the DAs of Ref.[19| and Refs.jy, [7| in 
FIGs. (|4l5[) . which are shown by the dashed lines and the 
dotted lines respectively. As for in different to its 
asymptotic behavior, it will provide great suppression in 
the end-point region. Such a behavior will be helpful to 
derive a reasonable power behavior for twist-3 DA [l7j . 
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TABLE V. Wavef unction parameters for tpp(x,U.±) with varying DA moments. 



(e 2 ) 


0.196 


0.248 


0.324 


<S> 


0.207 


0.262 


0.342 


0.207 


0.262 


0.342 


0.207 


0.262 


0.342 


A p (GeV- 2 ) 


93.4987 


92.4429 


90.9983 


89.7356 


88.7233 


87.2936 


84.7201 


83.7330 


82.3809 


Bp 


1.2835 


1.2926 


1.3056 


1.3177 


1.3261 


1.3385 


1.3679 


1.3756 


1.3868 


Cp 


1.3579 


1.4165 


1.5005 


1.3635 


1.4210 


1.5045 


1.3710 


1.4277 


1.5097 


MGeV) 


0.5849 


0.5882 


0.5930 


0.5950 


0.5984 


0.6033 


0.6095 


0.6131 


0.6181 



TABLE VI. Wavefunction parameters for ij)™(x,'k.±) with varying DA moments. 





0.077 


0.102 


0.127 




0.074 


0.094 


0.122 


0.074 


0.094 


0.122 


0.074 


0.094 


0.122 


A a {GeV~ 2 ) 


196.486 


194.804 


192.128 


187.336 


185.532 


182.777 


178.403 


176.532 


173.584 




-0.0284 


-0.0309 


-0.0349 


-0.0115 


-0.0144 


-0.0191 


0.0053 


0.0019 


-0.0032 


c a 


0.0922 


0.1108 


0.1364 


0.0885 


0.1069 


0.1322 


0.0845 


0.1026 


0.1276 


Pa{GeV) 


0.4107 


0.4123 


0.4147 


0.4169 


0.4187 


0.4213 


0.4234 


0.4253 


0.4283 




FIG. 4. Pion twist-3 DA The solid line, the dashed line 
and the dotted line are for our DA defined by Eq. (|12[l . DA of 
Ref.[li] and DA of Refs.0,0, respectively. 
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FIG. 5. Pion twist-3 DA <j>Z. The solid line, the dashed line 
and the dotted line are for our DA defined by Eq. (|13jl . DA of 
Ref.HU and DA of Refe.@,0|, respectively. 



It is noted that to study the right behavior of the DA it- 
self, the full form of the twist-3 DAs ([12")) and (TT3"f is more 



useful and more accurate than its truncated Gegenbauer 
form. 

B. reanalysis of the B — > tv transition form factor 

The B — ¥ 7T transition form factor provides a good plat- 
form to study the pion properties. Especially, within the 
QCD light-cone sum rules |33 |. one can concentrate on 
different twist structures of pion DA by properly choos- 
ing the correlator. Since by choosing the corrector with 
proper current, because of the cancelation among the dif- 
ferent 7-structures, different twist DAs will be remained 
in the final formulae. 

Three typical correlators are suggested in the litera- 
ture to calculate fg w (q 2 ). The first one is IL^(p,q) = 

i J d 4 xe^ x (n{ P )\T{q(x)^b(x),b(0)irnbJ 5 q(0)}\0), which 
is usually adopted in the literature. After doing the 
simplification, it can be found that all the twist-2 and 
twist-3 and twist-4 terms are remained and all of which 
can provide sizable contributions [3oT - l37j . The sec- 
ond one is 11+ (p, q) = if d 4 xe tq ' x ' {■n{p)\T{q{x)^ lJi {l + 

j5)b(x), b(0)imt)(l + 75)g(0)}|0), which is constructed 
by using the chiral current. After doing the sim- 
plification, it can be found that the twist-3 terms 
can rightly be canceled, so one only needs to con- 
sider the twist-2 and twist-4 DAs [H, HI]. The 
third one is Il+(p,q) = i J d A xe l ^ x {-K{p)\T{q{x)-i^{\ + 

75)6(12;), b(0)imt)(l — 75)g(0)}|0). After doing the simpli- 
fication, it can be found that the terms involving the 
twist-2 DA are exactly be canceled, so one only needs to 
consider the dominant twist-3 DAs [40|. 

In the present paper, our purpose is to test the proper- 
ties of the obtained twist-3 DAs in the last section, so we 
adopt the third correlator to do our analysis. Following 
the same procedure of Refs.[!3, 53 > we can obtain the 
light-cone sum rules for the /s w (9 2 ) U P to NLO, 
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jbM ) = —2-r e ^p 



~m B ~ 




M 2 _ 





Fo( (Z 2 ; M 2 , S o) + ^F 1 (g 2 ,M 2 , So ) 



(16) 



r 



where the leading order Fo(q , M , So) takes the form 



F (q 2 ,M 2 ,s ) =m b 
2JW 



/ttMtt Jo 



du 



■ exp 



to 2 — q 2 (l - u) 



A U 
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vdv / Da 



(ai + va 3 ) 



uM 2 
1 ex P 



m 2 - g 2 (l - ax - va 3 ) 



(ai + va 3 )M 2 



? 2 

- q l 



(ai + va 3 )M 2 



<h*(oi) K17) 



r 



and the NLO Fi(q 2 ,M 2 , sq) takes the form 



^(g 2 ^ 2 ,^) 



Trrrib 



ds e- s ' M ' 2 / [lm fl Tf (g 2 , s, u)0J(«) + Im s T[(q 2 ,s, u)^(u)] 



(18) 



where -Dai = e?aicZa2<ia:3<5(l — — a 2 — 013), mj, stands 

for the running b-quark mass under the MS scheme, A = 

2 _ 2 

™ b _ g g 2 , </>37r(ai) is the pionic three-particles twist-3 DA 
with the normalization parameter f3^ . The leading order 
result agree well with that of Ref.[40|. And it has been 
found that our present NLO result is different to that of 
Ref.[13] only by an overall factor '2', which is caused by 
the use of different correlators. So to short the paper, 
we do not present them here, and the interesting readers 
may turn to Ref.[33| for detailed formulae. Moreover, 
the three-particles twist-3 DA can be derived from the 
exact relation between the two-particle and three-particle 
twist-three DAs, which can be derived by the equations 
of motion 19]. It is noted that the contributions from 
<j>3Tr(c<i) to the present sum rules ([To]) is quite small, which 
is less than 0.5%, so we simply take its form as suggested 
in Ref.fl9j to do our calculation. 

We present our results for fg n in FIG.©, where the 
light-cone sum rule with certain extrapolation [39l |. and 
the quenched lattice QCD result [4l| and the unquenched 
lattice QCD result 42[ are also presented for comparison. 
The light-cone sum rule for B — > tt is valid up to momen- 
tum transfer q 2 ~ to 2 — 2mbA, typically at < q 2 < 14 
GeV 2 , and to be on safe side, we take the maximal al- 
lowed q 2 to be 12 GeV 2 . Here in doing the numerical 
calculation, the center value of f t r i s obtained by taking: 
f B = 169 MeV, s = 32.6 GeV 2 [39j and fh b {fh b ) = 4.20 
GeV, hi b — 5.279 GeV fH. The effective threshold for 
different correlators should be differ from each other so 
as to derive a reliable estimation [43[ . As for the present 
adopted correlator with chiral current, the value of sq 
should be slightly lower than that of the conventional 
corrector (i.e. the first correlator) so as to include the 



4p 
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3 - 

cr 

""K 2 

+ m 

II 1-5- 
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FIG. 6. /a 7r (9 2 )- The shaded band is our result with uncer- 
tainties caused by the input parameters, and as a comparison, 
the extrapolated light-cone sum rules of Ref.[39l] is shown by 
a dashed line. The quenched [4lJ and unquenched lat- 
tice QCD results are presented by rectangles and triangles 
respectively. 



unwanted scalar state with J p = + into the continuum 
contribution [Hj]. 

It is found that our present result agrees with those 
of Refs.jy, [jj, [l9| under the same parameter values. The 
shaded band is obtained by varying the pionic parameters 
within their reasonable regions. As for the Borel param- 
eter, its reasonable region is [14, 20] GeV 2 and we fix its 
value to be M 2 = 16 GeV 2 to do our calculation. The 



typical energy scale for B — > it is fi b = (m 



8 



2.2 GeV, and by running [/,„ from 1 GeV to (Mb, we ob- 
tain fi w ([i b ) = 1.93 ± 0.24 GeV. At q 2 = 0, we obtain 
/^(0) = 0.253 ± 0.031. It is found that the mam uncer- 
tainty is caused by fi^, and the twist-3 DAs models ([12")) 
and (|T3|) lead to small uncertainties. As for the NLO cor- 
rections, it will provide sizable contributions and its vaule 
will increase with the increment of q 2 , e.g. for q 2 = its 
contribution is ~ 6% and for q 2 — 10 GeV 2 its contri- 
bution changes to ~ 16%, so it should be included to 
provide a sound estimation. As a cross check, it is found 
that our present result of fg v (0) is consistent with those 
of Refs.[36j-|39| within reasonable parameter regions. 



V. SUMMARY 

We have studied the twist-3 DAs and <\>% within 
the QCD sum rules. And by adding all the uncertain- 
ties in quadrature, the uncertainties for the moments 
are about 20% — 35%. Further more, these moments 
can provide helpful constraints on the twist-3 DAs and 
hence the wavefunctions. The present constructed twist- 



3 wavefunctions shall be useful for the kx factorization 
approach [l3| . where the transverse momentum depen- 
dence both in the wavefunction and the hard scattering 
part should be treated on equal footing. A consistency 
check of our twist-3 DAs is done by studying the B — > it 
transition form factor f£ with the QCD light-cone sum 
rules up to NLO. 

Recently, Ref. [45] points out a non- negligible contribu- 
tion of higher-twist processes in large pt hadron produc- 
tion in hadronic collisions, where the hadron is produced 
directly in the hard subprocesses rather than by gluon 
or quark jet fragmentation. So a better understanding 
of higher twist wavefunctions or DAs shall be helpful for 
further studies on higher twist contributions. Especially, 
the forthcoming RHIC and LHC measurements will pro- 
vide further tests of the dynamics of large-p t hadron pro- 
duction beyond the leading twist. 
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